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Abstract 

We discuss the saturation of the parton density in heavy ion collisions at RHIC energies 
using a Pomeron approach. Our predictions for the particle density in ion-ion collisions at 
RHIC energies can be utilized as the background for the observation of possible quark-gluon 
plasma production. 
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For heavy ion collisions at high energy, saturation of the parton density is predicted both in high 
density QCD Jl|, || |J and in a Pomeron approach ||. The goal of this letter is to estimate the value 
of the inclusive cross section for the case of nucleus-nucleus interactions in a Pomeron approach, as 
this will be the background for any interesting signal of new physics such as the saturation of the 
gluon density |], 0, or the creation of a quark-gluon plasma at RHIC ||. The physical picture, 
behind our Pomeron approach, is the old fashioned parton model which was and can still be used 
as a guide in strong interaction physics at large distances ("soft" physics). In spite of the simple 
physics described by Pomeron exchanges, the Pomeron approach is technically rather complicated, 
and this problem has not yet been solved for hadron-hadron interactions. However, for inclusive 
production in ion-ion collisions there are two significant simplifications which enable one to develop 
a self consistent theoretical description in our Pomeron approach: 

1. For a heavy nucleus the vertex of the Pomeron-nucleus interaction gpA^t) — Aa gpN(b t = 
0) SU(&t) where gp^ipt = 0), is the vertex of the Pomeron - nucleon interaction at b t = 0, 
and SA(b t ) is the nucleus profile function normalized so that / d 2 b t SU^t) = 1. Due to Az 
enhancement of Pomeron-nucleus interaction we can develop a theoretical approach con- 
sidering GzpgpAipt = 0) ~ 1 while G\ P <C 1. Here, G 3P is the triple Pomeron vertex (see 
Fig. 1 ). In this approach one can neglect the loop Pomeron diagrams in comparison with the 
"tree" diagrams of Fig.l. 

2. Using the AGK-cutting rules ||, all Pomeron diagrams in which more then two Pomerons cross 
the rapidity level y c cancel in the inclusive cross section. Therefore, only Mueller diagrams, 
shown in Fig.l, survive. 

Using the approach, developed in Refs. |7|, |5], |9| ? we obtain the following closed expression for the 
diagrams shown in Fig.l ( see Fig.l for notation ) 



-- J d\ d\ a p P R Al (Y - y c , b t ) ■ R M (y c , b' t ) = (1) 

n P f J21. J2fJ gPAMgPAM) e AY 

p J 1 1 (gpAMiie^-y^-i) + i) {g PA M)i{e^-i) + 1) 

here 7 = G^p/A and G%p is the vertex of the triple Pomeron interaction, A = ap(0) — 1 where 
«p(0) is the Pomeron's intercept, a P is the particle emission vertex and Y = ln(s/lGeU 2 ) is the 
total rapidity interval for A\ — A2 collision. Ra{Y, bt) is the sum of "tree" diagrams and is given in 
il- 

Eq. ( |l|) predicts saturation of the density of produced particles which can be defined as 

A ATI \ da{A 1 +A 2 ) 

p( y c) = d -lM = *■ (2) 

dy c a tot (A 1 + A 2 ) 

where N{y c ) denotes the multiplicity of particles with rapidity y c . 

Eq.( |I|) leads to a constant density p(y c ) at high energies. To show this, we note that using a 
simplified Wood-Saxon type b t distribution for the density of a nucleus, namely, 

S a (h) = -^Q(R A - b t ) , (3) 
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Figure 1: The Mueller diagram for an inclusive A- A interaction process. 



the total cross section ( for A\ > A2 ) is 

a tot = 2rrR 2 Ai . (4) 
Integrating over b t and b' t in Eq.( |l|) we have 



R 2 

p{Vc) — ► |y»j/ c »i a P p TT^i > ( 5 ) 



where at high energy (see Fig.l) 



a P 



da(N+N) 
dye 



a tot (N + N) 

Eq. ( |) suggests that the particle density is independent of the energy and of the atomic number of 
the heaviest nucleus. It is instructive to compare this prediction with the density obtained in the 
Glauber-Gribov approach |I0| where 

P(Vc) — |„>i a p P A 2 Afe AY . (6) 

Eq.( ||) predicts a different behaviour for the particle density than does Eq.( H), with a clear 
dependence on energy and on A\. Eq. (|l]) also contains a very simple prediction for the particle density 
in the fragmentation regions. Assume that y c — > and y c — > Y correspond to the fragmentation 
regions for lightest and heaviest nucleus, respectively. From Eq. ([!]) it follows that 

p(Vc) | yc ->o = A 2 ; (7) 

P(yc)\y c ->Y = 4 4^- (8) 
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Figure 2: Our prediction for particle production. 



Fig. 2 illustrates our prediction for particle density. 

The physics that Eq. fl5|) is based on, is the same saturation of the parton density that has been 
discussed over the past decade in the framework of high density QCD |], |], ||. Namely, due to the 
emission of partons, the parton density increases at high energies but, simultaneously, the probability 
of the interaction of two partons increases as well. Since, two interacting partons can annihilate into 
one, such interactions diminish the number of partons in the parton cascade. Finally, competition 
between the emission and the annihilation leads to an equilibrium state with a definite value for the 
parton density. We call this phenomenon, saturation of the parton density. The high density QCD 
approach leads to the following A — A dependence of the process §], f|: 

< N Al A2 >= Q 2 S (s, Ai) -Al<N NN >, 

here A 2 < A 1 , Q 2 S (s, Ax) is the saturation scale, which increases with A, and the energy squared s, 
< N > is the density of the produced hadrons in the A\ — A 2 or the N — N interactions. Q 2 S (s, A) 
grows with A ( at least Q 2 S oc A~s |TI| , but it could even be proportional to As |T2J ), and therefore 
high density QCD predicts 

<n AiA2 >kaI-(a!+4), (9) 



(10) 
and 



>oc Al 

where A 2 < A\. Eq. (^) can be rewritten as: 

2 

< iV Al A2 >= Al <N NN >, 

where A 2 < A\. 

Eq. (0) leads to different predictions than those of both the Glauber approach (see Eq. 
the high density QCD approach ( see Eq. @ ). 
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To calculate the particle density at RHIC energies we need to specify how to calculate the total 
cross section. It turns out that even in our approach for nucleus-nucleus scattering, the problem of 
determining the total cross section is rather complicated, but it has been solved in Ref. 0. We use 
the following expression for o tot : 



crtot = 2 J d 2 b t ^1 - exp^-^^j j . (11) 

(12) 

Here fl(s, b t ) the opacity is given by 

Q(s,bt) = I d% F(Y,b' t ,\b t -V t \) , (13) 



and F [Y, b' t ,\b t — b' t \J is the amplitude of the nucleus-nucleus interaction |5|]. Since the RHIC energies 

are not very high ( Y = ln(s/lGeV 2 ) ~ 10 we can use a simple formula for F ( see Ref. || for 
details ): 

F (Yb b') g PAl (b t -V t )-g PA2 (V t )-e AY (u) 

[ ' " t} {gpM{h-bi t )+gPAM))l{e^ + l • 1 ] 

We assume the Wood-Saxon parametrization |T3| for profile function S A (b t ). 

S A (b t ) = (15) 

1 + e — ~ 

where R A is the nucleus radius R A = 1AA^ fm and h the surface thickness h m lfm. For a numerical 



estimate we used the parameters of Eq. fllq) given in Ref. |14 . 

It is even easier to calculate the particle density for central A-A collisions which has been measured 
at RICH. The observable can be written as a ratio of the inclusive cross section to the cross section 
at fixed bt = 0. 

r x J d%' t a p P R Al (Y - y b' t ) ■ R M (y c , b>) 

Pcentral {Vc) ~ —T. 7 nf^-nU \ l ib J 



2(1- 
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Q(s,bt=0) 
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d 2 h , 9paM) ■ 9paM) ■ e Ay (m 

1 (gpAMh (e A(y -^ - i) + 1) (9paM)i (e Ay < - 1) + 1) 1 ] 

since ft(s, b t = 0) ^> 1 for A-A scattering. 

We calculate the particle density, defined be Eq. (Q), for Au-Mo, Au-Ne and Ne-Mo scatter- 
ing at the RHIC energies, W = 56 GeV ( Y = In (s/1 GeV 2 ) = 8 ) and W = 156 GeV ( 
Y = In (s/1 GeV 2 ) = 10) . For proton-proton scattering we use the following parameters, which 
fit the experimental data quite well ( see Ref. f| for details ): 

7 = 1.19 GeV^ 1 g PN = 8AGeV- 1 A = 0.07 a p P = 2. 

Fig. 3 shows that we find a particle density which is 2 -=-3 times smaller than in the Glauber approach. 
The contrast is even more impressive when compared with the prediction of the gluon saturation 
approaches or/and for quark-gluon plasma production ( see Refs.Jl5| and for example ). The 



physics is very simple and is the same as for the gluon saturation approach: the interactions between 
partons diminish their number, but in contrast to gluon saturation the interaction only occurs when 
the partons have low transverse momenta. 



5 



W = 56 GeV 



W = 156 GeV 



dN/dy c 




dN/dy 



150 



100 



Glauber Formula 



Our Approach 



Au 




dN/dy c Glauber Formula 



Our Approach 




Au 



Mo 



Glauber Formula 



dN/dy c Q| au |j er Formula 






dN/dy, 




Glauber Formula 



Our Approach 



Au 



Au 




dN/dy c Glauber Formula 




dN/dy c 


250 




200 




150 




100 




50 





dN/dy c 


175 




150 




125 




100 




75 




50 




25 





Our Approach 



Au 



Mo 



Glauber Formula 



2 4 6 

Glauber Formula 



Ne 



Mo 



10 




10 




10 




10 



Figure 3: The density of produced hadrons for Au-Au, Au-Mo, Au-Ne and Mo-Ne in central A- A 
interactions at RHIC energies. 
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Comparison with the experimental datajnj[] is given in Fig.4. We plot three different theoretical 
curves in Fig.4: the upper one is the calculation of the hadron density for a central collision, calculated 
using Eq. ([16]) with b t = 0; the middle curve is calculated , assuming that centrality of the nucleus 
interaction was fixed with an accuracy b t < 2 fm Q and the lowest curve is the inclusive density of 
Eq. (Q). One can see that the agreement both in value and in the energy dependence is rather good, 
but it depends strongly on the accuracy with which the centrality of the experimental interaction 
was determined. However, it is very difficult to distinguish our prediction for gold-gold collisions 
from those based on the saturation approach at high density QCD, when the saturation scale is low 
(approximately 1 -j- 1.5 GeV ) as was suggested in Ref. ||16|| . For this case both approaches describe 



the same physics at RHIC energies and we need more data on different nuclei to differentiate between 
these two approaches. 

Our conclusions are simple. We show that the RHIC data on the rapidity multiplicity in the 
central A-A collision can be used to test the different models Jl|, ||, [| [0], IE, [16). We predict that at 



RHIC, densities will be rather small ( 400-500 particles per unit of rapidity for gold-gold interaction, 
see Fig.4 ) and these data support the saturation hypothesis, at least for sufficiently small transverse 
momenta for which our Pomeron approach is valid. Secondly, our calculations should be considered 
as a background for a signal of new physics such as saturation of the gluon density in high energy 
QCD or/and quark-gluon plasma production. Fig. 3 shows that the measurement of A-dependence 
will be very decisive in distinguishing between the saturation due to the Pomeron interaction and 
the high density QCD approach. 




Figure 4: Energy behaviour of the central hadron density. 
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